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ABSTRACT 



Seventeen-month time series (May 1989 - October 1990) of current, 
temperature and conductivity were obtained from 100, 350, and 500 m depth at 
site P2, located on the 800 m isobath off Point Sur, and one-year time series 
(May 1990 - May 1991) of the same variables at similar depths were obtained 
from site P3, approximately 25 km farther offshore on the 1800 m isobath. 

Results show that no net growth or decay of eddy potential energy (EPE) 
occurred at either mooring site during their respective deployment periods. At 
mooring P2, baroclinic instabilities within the water column were signaled by 
downgradient horizontal eddy heat fluxes that converted mean potential energy 
(MPE) to EPE at both 225 and 425 m. The dominant balance at 225 m was 
between mean flow advection (source) and upward eddy heat fluxes (EPE to 
eddy kinetic energy, EKE), with additional losses coming from downstream 
advection by the eddy flow. At 425 m, the dominant balance was between 
downgradient eddy heat fluxes (source) and downstream advection by eddy flow 
(sink). Unlike 225 m, vertical eddy heat fluxes at 425 m were a weak source 
(EKE to EPE) while mean flow advection was negligible. At P3, the net balance 
involved only downward eddy heat fluxes (source) and downstream advection by 
eddy flow (sink), as mean advection and MPE-EPE conversions were negligible. 

Analysis of energetic events within the time series of terms in the EPE 
equation did not reveal any canonical or common pattern which would explain 
the temporal means described above, but suggest the flow in this region is highly 
variable. In fact, during most events magnitudes of terms were anywhere from 



in 



•7 



10 to 200 times that of the associated temporal mean. Events at P2 involved both 
horizontal and vertical processes and had longer time scales (several days to 
weeks) compared to those at P3, which had much shorter time scales and 
appeared to involve vertical processes only. 
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I. INTRODUCTION 



A. THE CALIFORNIA CURRENT SYSTEM - A GENERAL 

OVERVIEW 

The large-scale atmospheric forcing in the eastern Pacific Ocean consists of 
the North Pacific (sub-tropical) high, the Aleutian low, and in summer the 
thermal low over the western United States. The North Pacific high is most 
intense during the summer months while the Aleutian low is most intense during 
the winter months. The high migrates annually from its maximum southern 
position at 28°N, 130°W in February to its maximum northern position at 38°N, 
150°W in August ( Huyer , 1983). The U.S. thermal low is centered near 35°N in 
summer and enhances the equatorward wind stress over the coastal waters off 
northern California ( Reinecker and Ehret , 1988). A region of positive wind 
stress curl exists near the coast throughout the year, is best developed from May 
to September, and has greater spatial variability during winter ( Nelson , 1977). 

This forcing creates the anticyclonic North Pacific gyre. The northern side 
of the gyre is comprised of the West Wind Drift and the North Pacific Current 
which flows easterly then splits near 45°N in winter, 50°N in summer ( Pickard 
and Emery, 1982) to form the poleward flowing Alaska Current and the 
equatorward flowing California Current (CC). Offshore the CC has been 
characterized as a surface current (0-300 m deep) carrying water equatorward 
throughout the year along the west coast of North America ( Lynn and Simpson, 
1987). Near 20°N it turns westward as part of the North Equatorial Current. 
The average speed of the CC off the coast of California is generally less than 25 
cm s* 1 ( Reid and Schwartzlose, 1962). 
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Within 150 km of the coast there is a fall- winter reversal of the surface flow 
known as the California Countercurrent (CCC, Simpson et al., 1986). This 
current has also been called the Inshore Countercurrent (IC, Lynn and Simpson, 
1987) off southern California, and is generally referred to as the Davidson 
Current (DC) north of Point Conception. We adopt the term Davidson Current 
(DC) to refer to this feature off Point Sur. The reversal of winds from 
northwesterly in summer to southeasterly in winter, which causes downwelling at 
the coast, seems to be the forcing mechanism of this poleward surface current 
( Huyer et al., 1989). 

The California Undercurrent (CUC) flows poleward throughout the year. 
While it has been observed off northern Baja California (Wooster and Jones, 
1970), central California (Chelton, 1984; Lynn and Simpson, 1987; Tisch et al., 
1992), Oregon ( Halpern et al., 1978), Washington ( Cannon et al., 1975; Reed 
and Halpern, 1976), and as far north as Vancouver Island, British Columbia 
(Reed and Halpern, 1976; Ikeda et al., 1984), the spatial continuity of this 
current along the west coast of North America has not been observed. It has its 
origin in the eastern equatorial Pacific and is centered primarily over the 
continental slope. Based upon available observations of the winds and the CUC 
between San Francisco and Baja California, Hickey (1979) concluded that the 
location, strength and core depth show considerable seasonal variability and can 
be related to the seasonal variability in wind stress and wind stress curl. Similar 
results were drawn from simple correlations between the current and variability 
within these fields (Lynn and Simpson, 1987) although many of the dynamical 
features have not been completely justified (Chelton, 1984). Tisch et al., (1992), 
using hydrographic data collected off Point Sur, California, found the location, 
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strength and core depth to be strongly related to specific wind events, both local 
and remote. While these observations have revealed that the core depth and 
speed can vary from location to location, little is known about what causes the 
observed variability and more importantly what actually drives the undercurrent 
(Huyer et al., 1989). 

B. SUBTIDAL VARIABILITY WITHIN THE CALIFORNIA 
CURRENT SYSTEM 

Collectively the CC, DC, and CUC comprise what is known as the California 
Current System (CCS). Much of our present understanding of the circulation 
within the CCS has come from analyzing hydrographic data, moored current 
meter measurements, buoy drift trajectories and, more recently, satellite imagery 
of the sea surface, which has been collected in conjunction with observational 
programs conducted within the CCS. A common goal of these programs has 
been to study the kinematics and dynamics within the CCS and eastern boundary 
current regions in general, and to examine the role played by these currents in 
larger gyre-scale circulations. 

Of particular interest in some of the more recent studies are: 1) the time 
variability of poleward flows and their role in gyre-scale processes; 2) the 
dynamical processes which govern the wind-driven circulation over the 
continental slope; and 3) the nature and structure of cold filaments found in 
eastern boundary current regions. In addition to the aforementioned 
observational programs, several numerical models of the CCS have been 
developed to further investigate possible mechanisms for producing the observed 
variability in eastern boundary currents. Using simple three-dimensional linear 
models forced by realistic winds from the region off central California, 
McCreary (1981) and McCreary et al ., (1987) have been successful at producing 
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poleward undercurrents. The poleward surface flow and undercurrents 
generated by these models depended upon the alongshore wind stress, and the 
existence of an alongshore pressure gradient. Models forced solely by curl of 
the wind stress did not produce poleward undercurrents, however, they did 
produce poleward surface currents nearshore. The primitive equation model of 
Batteen et al. (1989), forced by a band of steady alongshore wind, also generated 
an equatorward coastal jet and poleward undercurrent, which became unstable 
with time, and led to the production of eddies and jets with significant onshore 
and offshore directed flows. In a similar experiment, Batteen et al., (1989) also 
found that a variation of alongshore wind stress can play a role in determining 
the location of eddy generation regions. These studies, however, are just part of 
a vast amount of literature on the subject of eastern boundary regions, and more 
specifically, the CCS. For the sake of brevity, the following discussion will 
focus on descriptions of the variability observed within the CCS during some of 
the large-scale systematic programs, with particular attention to the time and 
space scales of these motions. 

Chelton (1984) and Lynn and Simpson (1987), using 23 years of 
hydrographic data collected in conjunction with the California Cooperative 
Fisheries Investigations (CalCOFI) program, examined the seasonal variability of 
alongshore geostrophic currents along the west coast of the United States. This 
large-scale sampling grid was initiated in 1949 and extends from Vancouver 
Island near the Canadian border to the southern tip of Baja California. These 
historical data reveal the CC within the upper 200 m flowing equatorward with 
its core located between 100 km and 200 km off the coast. The core location 
appears to be coincident with a transition zone that separates the coastal region 
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from the oceanic environment and is characterized by recurrent eddies and 
energetic meanders ( Lynn and Simpson, 1987). There are two maxima during 
the year, between February and March and July through September 
(climatologically the time of maximum equatorward wind stress) where the 
velocities exceed -9.0 cm s* 1 . The seasonal average geostrophic flow in the upper 
100 m relative to 500 m off Point Sur and Point Conception is equatorward flow 
from February to September and poleward flow from October to January 
0 Chelton , 1984). The CUC was present over the continental slope for most of 
the year, with the exception of early spring (March - May) and poleward flow 
extended to the surface in the wintertime (October - February). The surface 
flow throughout this region was found to lead the annual wind forcing, computed 
from the spatial averages of Nelson (1977), by about one month while the deeper 
poleward flow was found to lag the poleward barotropic pressure gradient by 
about two months ( Chelton , 1984). These time scales are in reasonable 
agreement with the modeling results of Philander and Yoon (1982). Using 200- 
day periodic alongshore wind stress, they found that the alongshore surface 
velocity led the alongshore pressure gradient and wind stress by approximately 
25 days, while the undercurrent lagged the alongshore pressure gradient by 
approximately 75 days. 

Wickham et al. (1987) examined current meter and hydrographic data 
collected off Cape San Martin, California between 1978 and 1980 and found that 
the undercurrent had a jet-like core with maximum velocities in excess of 15 cm 
s* 1 , and was confined to within 30 km of the coast in the upper 300 m. The 
currents exhibited a strong annual cycle, similar to the results of Chelton (1984). 
Tisch et al., (1992) found that the position of the core off Point Sur, California 
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varied between 12 and 42 km from shore, and between depths of 70 and 460 m, 
with speeds ranging from less than 5 cm s' 1 to 35 cm s* 1 , with maximum flow 
occurring in winter. 

While these observational studies began to identify the kinds of variability 
present within the CCS, and some of the basic time scales, the Coastal Ocean 
Dynamics Experiment (CODE) was the first study to address the issue of 
identifying and studying the dynamical processes of the wind-driven motion of 
water over the slope. This experiment was one of the first heavily instrumented 
programs of its kind conducted off the west coast of the United States. CODE 
was located off the northern California continental shelf between Point Reyes and 
Point Arena and occurred during the spring and summer of 1981 (CODE 1) and 
1982 (CODE 2) ( Beardsley and Lentz , 1987). Results from CODE 1 showed 
that the current field over the shelf was highly coherent in the vertical but 
exhibited significant horizontal mesoscale variability suggesting that the flow 
over the shelf and slope may be strongly influenced by eddy- like features, both 
offshore and within the coastal flow, local and remote wind forcing, and local 
topographic effects {Beardsley and Lentz, 1987). Therefore, in CODE 2, high 
vertical resolution was replaced with increased horizontal instrumentation to 
investigate this mesoscale variability. 

Strub et al. (1987) investigated the annual variations in atmospheric forcing, 
currents, water temperature, and sea level along the west coast of the United 
States between 35° N and 48° N. They found that monthly mean winds in 
fall/winter were poleward north of 35° N to 45° N for three to six months, and 
were weakly equatorward or zero south of 35° N. The monthly mean 
alongshore currents over the mid-shelf and shelf break for depths of 35 m or 



6 



deeper were poleward, and were associated with higher coastal sea levels and 
relatively warmer water temperatures. In spring/summer, they found the 
monthly mean winds to be equatorward for three (near 48°N) to six months 
(near 35°N). At these times sea levels were lower and the water temperature 
cooler. The monthly mean currents at 35 m over the shelf were equatorward 
from one to six months, being both longer in the north and over the shelf break 
than to the south or over mid-shelf ( Strub et al., 1987). The seasonal cycles of 
all variables showed a poleward propagation, with stations in the south leading 
those in the north by one to two months. Annual mean currents over the shelf 
break at 35° and 48° N were found to oppose the annual winds, while between 
35° N to 43° N, strongly fluctuating currents existed both in summer and winter 
regimes. 

Noble et al. (1987) found subtidal current fluctuations (33 hours < period < 
32 days) over the continental slope off northern California to be strongly 
polarized along local topography, and poleward in the mean. Alongslope 
currents on the upper slope were highly coherent in the vertical, with the first 
mode EOF (empirical orthogonal function) explaining approximately 15% of 
total variance. The alongslope variance was between 3 and 50 times greater than 
cross-slope variance. These observations also imply that the upper slope flow is 
mainly an extension of the poleward undercurrent observed over the California 
slope ( Chelton , 1984; Tisch et al., 1992), although as mentioned earlier, the 
continuity of the CUC has not yet been established. Mid-slope (instrument 
located on the 2200 m isobath) flow was generally weak, with no significant 
mean flow. Subtidal currents in the adjacent ocean basin were found to have no 
stable orientation, with a weak southerly mean flow. In this location, the easterly 
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variance was one to two times the northerly variance. The subtidal current 
fluctuations in each of these regions (upper slope, mid-slope, adjacent ocean 
basin) were generally incoherent with each other, and only weakly coherent with 
wind stress and sea level, whereas Winant et al. (1987) found that the subtidal 
currents over the shelf are, in general, coherent with sea level and wind stress. 
As a result, these observations indicate that shelf and slope flow may not be 
correlated with each other even when separated in the horizontal direction by as 
little as 25 km. 

Currents over the central California shelf and upper slope were also studied 
as part of the Central California Coastal Circulation Studies (CCCCS) program, 
conducted between Point Conception and San Francisco from February 1984 
through July 1985 ( Chelton et al., 1988). Consistent with the findings during 
CODE ( Winant et al., 1987), the alongshelf current fluctuations were found to be 
highly correlated with local winds (response within 0.5 day) and propagated 
poleward in rough agreement with the second-mode coastal trapped wave 
estimated by the model of Chapman (1987). In both the CCCCS and CODE 
regions, however, the correlation between currents and winds beyond the shelf 
break showed a marked decrease, and illustrate that currents separated by only 
10-15 km cross shore show less correlation than shelf currents separated by as 
much as 220 km alongshore (Noble et al., 1987; Chelton et al., 1988). 

During the Northern California Coastal Circulation Studies (NCCCS) 
program, conducted between San Francisco and the Califomia-Oregon border 
from March 1988 through October 1989 ( EG&G , 1988), the average annual 
flow was generally more poleward than the seasonal average flow during CODE 
(. Magnell , 1991). Strong poleward flow was observed over the shelf at nearly 
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every instrument location during late summer and was in opposition to 
equatorward wind stress. Very low frequency current fluctuations, O(months), 
were observed at all mooring locations and found to be poorly correlated along 
the coast, and with the local wind stress. Currents in the 2 to 20 day band were 
found to be correlated best with winds to the south, suggestive of remote forcing. 
This is consistent with the findings of Davis and Bogden (1989) who, based on 
CODE data, discuss a length scale of 500 km associated with remote wind 
forcing. Magnell (1991) also found that alongshelf pressure gradients, forced by 
alongshelf variations in wind stress (iy), were highly correlated with the local 
current accelerations. 

Strong poleward mean flow, in opposition to the equatorward wind stress, 
was observed during the CCCCS experiment and spread offshore to a distance of 
300 km. Similar occurrences were observed in July 1981 ( Chelton et al., 1988) 
and in July 1989 ( Tisch et al., 1992), where poleward flow extended to 80 - 150 
km offshore. These bursts of poleward flow were believed to have been the 
result of large-scale wind relaxations that occurred previously over the central 
and southern California regions. The coastal response to a relaxation from 
upwelling were examined by Huyer and Kosro (1987) and Send et al. (1987) 
during CODE. While the appearance of poleward surface flow and associated 
warm water, in response to a cessation of equatorward upwelling favorable 
winds, is similar to these larger poleward bursts, the poleward flow observed in 
the CODE region appeared to be trapped to a narrower region over the 
continental shelf, and had a time scale of approximately one week (compared to 
nearly six months in 1981 and 1984) ( Chelton et al., 1988). 
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Utilizing data from the Ocean Prediction Through Observation, Modeling, 
and Analysis (OPTOMA) program, Reinecker et al. (1988) found strong vertical 
coherence at all mooring locations on the continental rise in the upper 600 m of 
water column, with the velocity components being in phase for periods longer 
than about 10 days (over vertical separations of up to 500 m). Motions in the 
upper ocean were not highly coherent with those in the lower ocean, and band- 
averaged coherences of u, v, and T as a function of vertical separation showed a 
decrease in coherence as vertical separation increased for periods longer than 10 
days. Temperature was found to be less coherent at smaller separations than u 
and v. At 350 m depth, the highest coherence was found between the out-of- 
phase v components at periods of 45-90 days. At periods of 9 - 13 days, they 
observed a surface and bottom intensification of the nearshore kinetic energy, 
which they found to be partially consistent with the presence of topographic 
Rossby waves. 

Adding to the mesoscale variability within the CCS are filaments and jets, 
which extend offshore from the shelf break to several hundred kilometers ( Brink 
and Cowles , 1991). Filaments are commonly observed features in satellite sea 
surface temperature and ocean color (CZCS) imagery of the central California 
coastal waters from Cape Blanco to Point Conception ( Bernstein et al. 1977; 
Breaker and Gilliland, 1981; Ramp et al., 1991a). It has been the goal of the 
Coastal Transition Zone Experiment (CTZ, conducted from 1986 through 1988 
offshore of the CODE region), to study, among other things, the physical 
structure and characteristics of these cold filaments, with an emphasis on 
furthering the present understanding of the kinematics and dynamics involved. 
Several theories have been proposed for their physical cause; 1) baroclinic 
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instability of the flow within the CCS ( Ikeda and Emery , 1984); 2) irregular 
coastline geometry ( Peffley and O'Brien, 1976; Crepon et al., 1984); 3) variable 
bottom topography ( Preller and O'Brien, 1980); 4) alongshore variation of wind 
stress ( Batteen et al., 1989): and 5) interaction with an offshore eddy field 
( Reinecker et al., 1988; Reinecker and Mooers, 1989). However, which cause(s) 
will generate these features in a particular geographical location is still an area of 
active research. 

Analyzing buoy drifter tracks from the CTZ experiment. Brink et al. (1991) 
found that in the summer and fall, the CC may be characterized as a meandering 
coherent jet, which on average flows southward to at least 30° N. Between 33° N 
and 39° N, core velocities of 50 cm s' 1 were typical, and current meanders had 
alongshore wavelengths of 0(300 km) and cross-shore amplitudes of 0(100-200 
km). These observations are consistent with those of Bernstein et al. (1977), 
who observed the CC to be a meandering jet, with wavelengths of 300 to 500 
km. Brink et al. (1991) also found that this meandering can lead to large eddy 
kinetic energies and eddy diffusivities, especially north of 36 0 N. In general, the 
results from the CTZ program reveal that the core of the CC in summertime is 
unstable, leading to the formation of a meandering jet surrounded by persistent 
eddies {Brink and Cowles, 1991). Huyer et al. (1991) conducted repeated 
mesoscale surveys throughout the CTZ region and found that a meandering 
baroclinic equatorward jet was common to all. The core velocity exceeded 50 - 
70 cm s' 1 , and may be coincident with the core of the CC. They also observed a 
poleward flowing undercurrent adjacent to the continental slope at depths of 150 
- 250 m, with core velocities up to 20 cm s* 1 . 
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From these studies we have begun to assimilate a vast amount of information 
on the kinematical and dynamical aspects of the CCS, and eastern boundary 
regions in general, and have demonstrated that the CCS is not a simple, but 
rather a complex and highly variable flow regime, consisting of eddies, 
filaments, and energetic baroclinic jets. In light of this it is surprising that the 
study of energetic processes within eastern boundary regions has not gained 
more attention. If the poleward undercurrent, found in eastern boundary 
regions, is unstable, it is likely to play an important role in the generation of 
these cross-shore jets and eddies ( Mooers , 1989). Therefore, the study of 
energetics within eastern boundary regions is crucial to determine if this transfer 
mechanism is substantial, and to identify what energy conversion processes and 
balances actually exist. 

C. ENERGETICS OF OCEANIC GYRES, INCLUDING BOUNDARY 
CURRENTS 

The study of energetics of the world oceans has primarily focused on 
determining the effects of eddies and time-dependent phenomena on the mean 
circulation of oceanic gyres, including boundary regions. As discussed earlier, 
this eddy variability can act as a signal for such oceanic processes as instabilities 
associated with strong currents, <and as a mechanism for the transport of heat and 
momentum. The distribution of eddy kinetic energy throughout the world ocean 
has also provided a means for validating numerical ocean models (. Holland and 
Schmitz , 1985; Schmitz and Holland, 1982, 1986; Hall, 1991). 

In recent years, a great deal of effort has been put forth in studying the 
interactions of eddies and western boundary currents (WBC), largely because of 
the role that WBCs play in determining the nature of basin-scale flow (Dewar 
and Bane, 1985), and because the strongest oceanic eddy variability is associated 
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with WBCs, and their extensions (Hall, 1991). In-depth analyses of the 
energetics of eddy-mean flow interactions have been made for the Florida 
Current (Schmitz and Niiler, 1969; Brooks and Niiler, 1977), the Gulf Stream 
(Webster, 1961, 1965; Watts and Johns, 1982; Dewar and Bane, 1985, 1989a, b; 
Hall, 1986a; Rossby, 1987), and the Kuroshio Extension (Nishida and White, 
1982; Hall, 1991). The following discussion highlights some of the methods and 
key results obtained from some of these studies. 

Using data from multiple current meter moorings, located within the Gulf 
Stream (GS), Dewar and Bane (1985, 1989a, b) have computed energy budgets 
for the mean and eddy flow fields. This was accomplished by using the ensemble 
averaged horizontal momentum equations, multiplied (vector multiplication) by 
the ensemble averaged velocity. They found that a release of mean kinetic 
energy by the eddies constituted the dominant form of energy conversion within 
the GS off Charleston, South Carolina (Dewar and Bane, 1985), and that eddy 
pressure work may be important in the fluctuating energy budget. Their 
calculations also revealed that the mean kinetic energy flux within the GS in the 
South Atlantic Bight increases downstream, and the eddies tended to decelerate 
the mean flow. In order to obtain a balance in the mean flow equations, they 
concluded that the GS in the Bight must be releasing mean potential energy by 
flowing down a mean pressure gradient. Similar results were obtained in the GS 
farther to the north (Dewar and Bane, 1989a), where flow below 880 m depth 
was down a mean pressure gradient; however, above this depth, they concluded 
that a conversion of mean kinetic energy to mean potential energy occurred via a 
flow up a mean pressure gradient. 
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Hall (1986a) examined the energetics of the GS at 68° W, using eddy kinetic 
and potential energy equations, derived from the horizontal momentum and heat 
equations, respectively. Results indicated a net conversion of mean to eddy 
kinetic and potential energies, with energy exchanges being dominated by 
motions in the upper 1000 m of the water column. The well-defined structure of 
horizontal and vertical velocity, and temperature (Hall, 1986b) made it possible 
for this study to deduce the time-averaged structure from data from a single 
current meter mooring (Hall, 1986a). 

A similar study was conducted by Hall (1991) in the Kuroshio Extension, 
using data from a single current meter mooring. Once again, horizontal 
derivatives within the momentum equations could be estimated based upon the 
cross-stream structure of the Kuroshio Extension defined by Hall (1989a). In 
this study baroclinic and barotropic energy conversions within the current were 
examined using both geographic and "stream" coordinate systems. Cross-stream 
position was quantified in terms of the measured temperature at 350 dbar, and an 
alongstream flow direction defined by the measured current shear. Based upon 
this information the time series of velocity and temperature were rotated into a 
"stream" coordinate system, whose axes were attached to the meandering 
current. Using stream coordinates, Hall (1991) found significant conversions of 
mean to eddy potential energy on the anticyclonic side of the current, and 
smaller conversions of eddy to mean energy over the cold portion. In the 
geographic coordinate system, the Reynolds stresses were much stronger than for 
the stream coordinate system and resulted in relatively large apparent eddy to 
mean kinetic energy conversions. 
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While these and other studies have examined the nature of energetic 
conversions within western boundary regions from different perspectives, there 
is very little information about conversion processes within eastern boundary 
regions. As the results from the observational and modeling studies within the 
CCS have demonstrated, a great deal of mesoscale variability exists, and one 
could expect intense eddy-mean flow interactions to occur. To date no 
comparable studies of energetics over the continental slope within eastern 
boundary currents have been conducted, nor have there been any based upon 
direct observations within the CCS. Recently, however, a wind-driven eddy- 
resolving limited-area quasi-geostrophic (QG) model has been used by Auad et 
al., (1991) to examine the circulation and energetics within the CCS. The model 
includes quasi-realistic bottom topography, true coastline, and is embedded 
within a coarser model covering most of the North Pacific Ocean. QG theory is 
an approximation to the primitive equations when the Rossby number is small 
and the scaled solution remains 0(1) ( Walstad et al., 1991). Based upon 
comparisons between observations and modeling results within the Gulf Stream 
and Kuroshio Extension ( Schmitz and Holland, 1982, 1986), QG models are 
capable of accurately reproducing the observed eddy energy fields and their 
associated time scales ( Auad et al., 1991). 

The model simulated the main components of the CCS, namely the CC, CUC, 
DC, and a feature known as the Southern California Eddy, within the Southern 
California Bight. The simulated CC consisted of active eddy-mean flow 
interactions, and was an effective source of both first-mode baroclinic annual 
Rossby waves between 25° and 33° N, which propagated westward through the 
model domain, and long period ( > 200 days) waves north of 34° N, believed to 



15 



result from baroclinic instability (BCI). This eddy-mean flow interaction is 
consistent with earlier observations, which have characterized the CC as a 
meandering feature rich in mesoscale variability ( Lynn and Simpson, 1987; 
Brink and Cowles, 1991). Mesoscale variability within the model was 
characterized by periods of 0(100 days), and wavelengths of 0(200 km). These 
length scales are somewhat smaller than the observations of Bernstein et al., 
(1977) and Brink et al, (1991), but are consistent with the idea that the CC can be 
thought of as a meandering jet surrounded by eddies. 

Near-surface circulation was found to be energetically dominated by the 
mean flow, whereas the deeper (below ~ 3500 m) circulation was dominated by 
the eddy field. Energy transmitted by the wind to the mean flow in the upper 
500 m was both fluxed out of the open boundaries as a result of the P effect and 
also transformed into available potential energy (APE). Baroclinic processes 
were then responsible for converting this APE into a highly energetic eddy field 
(Auad et al., 1991). Below 500 m depth, energy appeared to be radiated 
vertically through action of the eddies themselves. In the core of the simulated 
CC, they found that downgradient eddy fluxes of temperature were of greater 
importance in the production of eddy kinetic energy than was the advection of 
available eddy PE. 

Walstad et al., (1991) recently investigated the dynamics of the Coastal 
Transition Zone off northern California, through the use of assimilation 
modeling techniques. In their study, a regional baroclinic QG model was driven 
by initial and boundary conditions derived from objective analysis of 
hydrographic and acoustic Doppler current profiler (ADCP) data, collected in 
late May and early June 1987. During this period, the CTZ consisted of a 
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meandering jet in the southwest portion, inflow from a jet in the northwest 
portion, a cyclonic eddy in the northeast portion and a poleward current in the 
southeast portion. Pierce et al. (1991) found that local Rossby numbers 
calculated across this meandering jet approached a maximum value of nearly 0.2, 
and concluded that the QG approximation still remains valid in this region. 

In their analysis, Walstad et al., (1991) found that the divergence of 
horizontal advection of kinetic energy was a mechanism forcing the jet meander, 
with pressure work redistributing energy both horizontally and vertically, while 
the total pressure work divergence nearly balanced the total time rate of change 
of kinetic energy following a parcel (DKE/Dt). They also found that the initial 
meander pattern which developed in the jet was primarily the result of 
barotropic instability (BTI), while buoyancy (the PE to KE conversion term) and 
vertical pressure work divergence were active on wavelengths greater than 200 
km. However, below the main thermocline in the vicinity of the meander, both 
baroclinic and barotropic mechanisms were found to be important. Based upon 
these results, they concluded that a large mixed instability meander was 
responsible for the propagation of the jet away from the coast and out of the 
region. 

The linear stability of a CT^ jet was examined by Pierce et al. (1991) using a 
six-layer QG model with observed velocity profiles, based upon the objective 
analysis of Walstad et al. (1991), serving as the basic state for the model. Their 
results indicated that barotropic instability within the cross-jet component of the 
divergence in advection of kinetic energy led to a net transfer of mean kinetic to 
perturbation kinetic energy. When perturbation wavelengths were less than 90 
km, the resulting instabilities were almost entirely barotropic in nature. For 
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wavelengths between 110 - 190 km, eddy potential energy began feeding back to 
the mean; however, the barotropic conversion process was still able to maintain 
the instability. At these wavelengths, there was also a strong transfer of eddy 
kinetic energy downward through the model layers. When the perturbation 
wavelength was equal to 260 km, both mean potential and mean kinetic energy 
were feeding the perturbation at all levels, signifying that both barotropic and 
baroclinic processes were contributing to the growth of the disturbance. The 
kinetic energy transformation was the dominant process in the upper 500 m 
depth, while below this depth the potential energy conversion exceeded the 
kinetic energy conversion ( Pierce et al., 1991). The 260 km wavelength was the 
fastest growing mode, with the ratio of BCI/BTI being nearly 0.9, while for 
wavelengths greater than 260 km, they found that the instabilities were more 
baroclinic in nature ( Pierce et al., 1991). 

D. SCIENTIFIC OBJECTIVE OF THE DISSERTATION 

The purpose of this dissertation will be to examine the energetics of the low 
frequency variability observed within the CCS using current meter data collected 
off Point Sur, California. To date, no detailed study of energetic interactions 
within eastern boundary currents, based upon direct observations, has been made 
over a continental slope region. A detailed study of the energetic interactions 
between mean and eddy fields will not only provide a means for comparison with 
the modeling results of Auad et al. (1991), Walstad et al. (1991), and Pierce et 
al. (1991), discussed above, but will hopefully provide us with a better 
understanding of the California Current System and eastern boundary currents in 
general. 
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While the energy studies of Hall (1986a, 1991) were based upon data from 
single current meter moorings, detailed descriptions of the cross-stream 
structure of the Gulf Stream and Kuroshio Currents provided the necessary 
information to evaluate several terms in the eddy kinetic energy (EKE) equation, 
as well as terms in the eddy potential energy (EPE) equation. The EKE 
equation, obtained by multiplying the u and v momentum equations by u’ and v’, 
respectively, adding them together and taking the time average, requires 
information on horizontal gradients of velocity. The data set used in this study 
does not contain enough information for the computation of these horizontal 
gradients. Furthermore a well-defined cross-stream structure for the CCS off 
Point Sur, comparable to those used by Hall (1986a, 1991), does not exist. 
Therefore, the method chosen to examine energetic conversions in this study 
must rely solely on the data set available, and will be limited to the use of the 
eddy potential energy equation. 

To accomplish this task, I have chosen an approach similar to that of Niiler 
and Hall (1988) and Hall (1991), in which they examined low-frequency eddy 
variability within the eastern North Pacific Subtropical gyre, and the Kuroshio 
Extension, respectively, using data from a single mooring. Niiler and Hall 
(1988) examined the energy conversions between the mean flow, very low 
frequency motion, and eddy band through use of the eddy potential energy 
equations. For their data set (a three year time series) the low frequency 
(periods greater than 200 days) and eddy bands (periods from 40 to 200 days) 
were determined from the spectral characteristics of the data. All computations 
in these two studies were performed using data which were synthesized at mid- 
depth levels between the instrumented depths. In creating these mid-depth time 
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series, Niiler and Hall (1988) and Hall (1991) were also able to obtain time series 
of the horizontal gradients of temperature (3T/3x, 3T/3y), which appear in the 
EPE equation, through the use of the thermal wind relations. Using the observed 
data from above and below, 3T/3x and 3T/3y were calculated at mid-depths 
using the velocity shear at the midpoint. 

In proceeding with this method, it has been necessary to synthesize time 
series of all variables (u, v, w, T) and VhT at the mid-depth levels between the 
actual instrumented depths. A complete description of the data sets and data 
processing methods, which includes the technique used to create each of these 
time series, is contained in Chapter II, with the exception of w, which is found in 
Chapter IV. General characteristics and details on the spectral analysis of these 
data sets can be found in Chapter III. The plan of the dissertation is to examine 
the energetic conversions through use of the eddy potential energy equation. The 
development of this conceptual model has been described in Chapter IV, with 
specific details being deferred to Appendix D. Time series of each term in the 
eddy potential energy equation have been computed using this model and allow 
for examination of energetics within individual events as well as in the mean. 
From this we may learn whether the mean values are the result of a few specific 
events or are actually representative of the mean state off Point Sur. Periods 
where locally intense bursts (growth/decay) occur within these terms will also be 
examined to determine how eddy potential energy is growing or decaying and 
what the dominant source and sink terms are. 
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II. DATA AND METHODS 



A. DATA COLLECTION 
1. Current Meter Data 

Vector-averaged current velocities, current directions, temperature, 
conductivity, and pressure were collected at two mooring sites, denoted P2 and 
P3, off Point Sur, California (Table 1). Mooring P2, located on the 800 m 
isobath over the continental slope, was equipped with three Aanderaa vector 
averaging current meters (Model RCM8) placed nominally at 100, 350, and 500 
m depths. Mooring P3, located near the 1800 m isobath, was equipped with four 
current meters placed nominally at 100, 350, 500, and 1000 m depths (Figure 1). 
Data were recorded at 30 minute intervals using solid state Data Storage Units 
(DSU), each of which is capable of storing 65530 10-bit words. At this interval 
each DSU could hold up to seven months of data, thus requiring data collection 
to occur over several deployments between May 1989 and May 1991 (Table 2, 
Figure 2). 

Current velocities were measured through the use of a shrouded paddle 
wheel and electronic counter assembly located at the top of the recording unit. 
Velocity, in units of cm s' 1 , was determined from the number of revolutions 
made by the paddle wheel during each sampling interval. The sensor had an 
accuracy of + 1.0 cm s _1 over its range of 2.0 - 250 cm s _1 , while below a 
velocity threshold of 2.0 cm s _1 the paddle wheel would stall ( Aanderaa , 1990). 
Both pressure sensors and the mooring dynamics program in use at the Naval 
Postgraduate School (NPS) indicated a maximum mooring tilt of 11.7° 
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Table 1. POINT SUR MOORING LOCATIONS 



Mooring 


Geographical Position 


Instrumented Levels 


Depth 


P2 


36 20.0' N 122° 10.2' W 


100, 350, 500 m 


800 m 


P3 


36 20.0’ N 122° 21.6' W 


100, 350, 500, 1000 m 


1800 m 



Table 2. CURRENT METER DEPLOYMENT DATES 



Mooring 


No. 


Deployment 


Recovery 


P2 


1 


11 May 1989 


24 August 1989 


P2 


2 


25 August 1989 


14 December 1989 


P2&P3 


3 


15 December 1989 


12 May 1990 


P2 


4 


14 May 1990 


9 October 1990 


P3 


4 


13 May 1990 


10 October 1990 


P2 


5 


10 October 1990 


12 May 1991 


P3 


5 


11 October 1990 


12 May 1991 



( Sielbeck , 1991), which is just below the acceptable level for proper compass 
operation; however, most of the time the mooring tilt was below this value. 

Current direction was measured by a magnetic compass/potentiometer 
assembly located below the storage unit. The compass was oil damped and 
required an average of 3 to 5 seconds to resolve directional changes of 90° or 
more. It could function properly up to a maximum tilt angle of 12° from the 
vertical and had a resolution of 0.35°. The accuracy of the compass varied 
between ± 5.0° for current velocities of 5 - 100 cm s* 1 , and ± 7.5° for velocities 
of 2.5 - 5 cm S ' 1 and 100 - 200 cm s' 1 . Combined magnetic deviation and 
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variation for each unit was determined on a surveyed test bench at NPS prior to 
each deployment and was accounted for in the postprocessing software ( Sielbeck , 
1991). 
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Figure 1. Current Meter Configuration at Moorings P2 and 

P3 (Adapted from: Sielbeck, 1991) 



2. Hydrographic Data 

During the deployments of moorings P2 and P3, oceanographic 
research cruises were being conducted as part of the Point Sur Transect (POST) 
and California Undercurrent (CUC) Programs conducted by the Department of 
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Figure 2. Current Meter Deployments at Moorings P2 and P3: 
The solid (blank) bars represent areas where good (no) data 
was obtained. The two shaded areas (P2-500m and P3-350m) 
represent areas where only u,v data was bad. Horizontal 
arrows denote individual deployments and the small vertical 
lines represent the times of hydrographic cruises listed in 
Table 3. 



Oceanography at NPS. The POST, established in 1987, has been occupied 6-8 
times yearly to examine long-term variability in the CCS off Point Sur, 
California. Data were collected at stations along the transect from the surface to 
within 50 - 150 m of the bottom using Neil Brown Instrument Systems (NBIS) 
Mark IIIB CTD's. For a complete description of the POST and CUC programs, 
along with detailed information on CTD calibrations, the reader is referred to 
Tisch (1990). Hydrographic data from stations located near moorings P2 and 
P3, collected during 14 cruises (Table 3), have been used to identify a 
relationship between temperature and pressure at each mooring and to compute 
density, p, and an equivalent compressibility coefficient, a, ( Dewar and Bane, 
1985; Niiler and Hall, 1988) at each instrument location. This information is 
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required for the calculation of terms found in the heat and energy equations and 
will be discussed in detail later in the chapter. 

B. DATA PROCESSING 

Upon recovery, raw data from the DSUs were downloaded to a computer 
and converted to standard scientific units using the sensor calibration equations. 
Data were then visually inspected for outliers or periods suspect of instrument 
failures or malfunctions, and were edited or deleted manually when necessary. 
The data were then filtered prior to further analyses using methods originally 
developed at Oregon State University ( Denbo et al., 1984). A standard right- 
hand coordinate system is used where x, y, and z are positive east, north, and up, 
respectively, and pressure, P, is positive downward. 

1. Filtering 

The postprocessed records were initially filtered with a Cosine-Lanczos 
filter utilizing a centered 25 point data window. This filter had a half power 
period of 2.9 hours (8.4 cpd) and completely removed signals whose periods 
were shorter than 2.0 hours. The 30 minute records were then interpolated to 
60 minute intervals using Lagrangian polynomials, which allow for interpolation 
of unevenly spaced data points to specified intervals ( Gerald and Wheatley , 1989) 
with each point falling on an even hour. At this stage any gaps in the hourly data 
sets were identified and filled using procedures described in the next section. 

2. Data Gaps and Gap Filling Procedures 

The presence of gaps in a time series has a negative effect on the 
proposed statistical analyses, especially when dealing with frequency-domain 
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Table 3. HYDROGRAPHIC CRUISES OCCURRING DURING 
CURRENT METER DEPLOYMENTS: Cruise 
designations beginning with ST were in conjunction with the 
POST program, while CU pertains to the CUC program. 
POST stations 3-6 and CUC stations 5, 6, 10, and 1 1 were 
used for this study. 



Cruise Designation 


Date of CTD casts 


STMAY89 


5/5/89 


CUJUL89 


7/29/89 


CUSEP89 


9/26/89 


CUNOV89 


11/16/89 


CUJAN90 


1/18/90 


CUMAR90 


3/2/90 


STMAY90 


5/2/90 


CUJUN90 


6/16-17/90 


CUAUG90 


8/13-14/90 


CUOCT90 


10/25-26/90 


CUDEC90 


12/13/90 


CUFEB91 


2/21-22/91 


CUAPR91 


4/5-6/91 


STMAY91 


5/8/91 



calculations. Therefore to obtain the longest possible segments, thus ensuring a 
continuous time line, gap filling techniques were required. Care must be taken 
to select a procedure, though ad hoc, which will introduce a minimum amount of 
error into subsequent analyses. There were two kinds of gaps in the data: short 



26 



gaps, usually between deployments, which could be filled using simple 
techniques, and gaps greater than 1 day which required more sophisticated 
methods. A period of 16 to 18 hours was common between recovery and re- 
deployment of instrumentation (Table 4). With the exception of the pressure 
time series, only those gaps in the hourly data sets where the gap length/record 
length ratio was 0.1 or less were filled. 

During the third deployment at mooring P2 (Table 2), the pressure 
sensor at the 350 m depth failed at deployment. To create a time series of 
pressure at this level, a value of 250 m, the amount of wire between instruments, 
was added to the pressure of the 100 m sensor. A similar procedure was used to 
correct for erroneous pressure values at the 500 and 1000 m instruments during 
the first deployment at P3. The time series at 500 m showed considerable noise, 
and the variance (a 2 = 81.86 dbar 2 ) was an order of magnitude larger than the 
upper level instruments (a 2 < 2.0 dbar 2 ). Such variability at depth is highly 
unlikely given the nature of flow throughout the region and is an indication that 
the sensor failed. The pressure at 1000 m showed less variability (a 2 < 9 dbar 2 ) 
than at 500 m; however it contained a drift over time. This drift was not present 
in data from the three instruments lying above it in the water column, indicating 
that the mooring's anchor did pot drag along the bottom and that this data was 
also in error. Here the time series from the 350 m instrument was used to create 
new time series at these levels by adding 150 m and 650 m, respectively, which 
once again was the amount of wire between instruments. This method does not 
account for the fact that the upper level instruments may be "blown over" more 
than the lower instruments, as is usually the case in stronger flow regimes such 
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Table 4. TEMPORAL EXTENT OF GAPS AT MOORINGS P2 
AND P3: The asterisks indicates that only u, v data was 
missing from the current meter record. 



Instrument 


Starting Date 


Ending Date 


Duration 


P2 - 100 m 


3 Jun 89/ 1200* 


4 Jun 89 / 0500* 


18 hrs* 




24 Aug 89/ 1800 


25 Aug 89 / 1000 


17 hrs 




14 Dec 89/0900 


15 Dec 89/0300 


19 hrs 




28 Apr 90/ 1600 


13 May 90/2300 


~15 days 


P2 - 350 m 


24 Aug 89 / 1700 


25 Aug 89/1000 


18 hrs 




14 Dec 89 / 1000 


15 Dec 89/0300 


18 hrs 




1 May 90/ 1500 


13 May 90/2300 


~13 days 




9 Oct 90 / 1500 


10 Oct 90/0800 


18 hrs 


P2 - 500 m 


23 Jul 89 / 0600* 


27 Jul 89 / 1400* 


105 hrs* 




24 Aug 89/ 1800 


25 Aug 89/1000 


17 hrs 




14 Dec 89 / 1000 


15 Dec 89/0300 


18 hrs 




12 May 90/2100 


13 May 90/2200 


26 hrs 


P3 - 100 m 


24 Apr 90 / 2000 


13 May 90/0400 


~20 days 




10 Oct 90/0900 


11 Oct 90/ 0100 


17 hrs 


P3 - 350 m 


16 Apr 90/2300* 


7 May 90 / 0100* 


~20 days* 




7 May 90 / 0200 


13 May 90/0400 


~7 days 




8 Oct 90/1100 


11 Oct 90/0100 


63 hrs 


P3 - 500 m 


7 Apr 90 / 0900 


11 Oct 90/0100 


~188 days 


P3 - 1000 m 


10 Apr 90/ 1100 


13 May 90/0400 


~34 days 




1 Jul 90 / 0900 


11 Oct 90/0100 


"102 days 



as the Gulf Stream and Kuroshio currents, where stronger currents (in excess of 
1.0 m s' 1 ) can cause vertical excursions of several hundred meters (Hall, 1989b). 
However, the small variability shown in the time series of pressure (Table 5) 
show that blowover was very small (less than 2 m) and would indicate that this 
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procedure can be used in this region without introducing significant error into 
subsequent calculations. 

Missing hourly data points of velocity (u, v) and temperature for the 
periods between deployments (gaps less than 18 hours) were interpolated by use 
of a cubic spline, while for pressure, which showed considerably less variability 
during deployments, simple linear interpolation was used. In these cases, the gap 
length/record length ratio was considerably smaller than 0.1, and the gap length 
itself was less than the half power period of the low pass filter. Since the focus 
of this study is on subtidal variability (periods longer than two days), it is 
expected that this technique will introduce minimal error into further analyses. 

To make the time series from mooring P2 complete for examining 
spectral properties and studying energetic conversions over a 17-month period, 
three additional gaps, one occurring at each of the instrumented depths and all 
less than two weeks in duration (Table 4), had to be filled using a more rigorous 
technique. Briefly, this technique consists of: l) initially filling the gaps with 
linearly interpolated values; 2) detrending the time series using least squares 
techniques; 3) estimating the sample periodogram; 4) removing any dominant 
cycles within the detrended data set (here being the diurnal and semi-diurnal 
signals): 5) fitting a first order autoregressive model, AR(1), to the detrended 
and deseasonalized data; 6) forecasting and backcasting from each end of the gap 
to interpolate the values across the gap; and 7) adding the estimated trend and 
cycle back to the interpolated times series. A similar procedure was also applied 
to the 63 hour gap in the 350 m record at mooring P3 (Table 4) so that a 
continuous record, one year in length, could be obtained. The specific details of 
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this technique, developed by P.A. Lewis and B.K. Ray of NPS, are contained in 
Appendix B. 

The gap-filled hourly data sets were filtered using a second Cosine- 
Lanczos filter, having a half power period of 46.59 hours (0.515 cpd) designed 
to remove diurnal, inertial, and shorter period energy within the signal. The 
filter utilized a centered 121 point data window to produce a single filtered data 
point. The filtered data were then decimated to six hourly values for subsequent 
analyses. Plots of current velocity vectors and temperature for moorings P2 and 
P3 using this filter are contained in Appendix A. Basic statistics of these data 
sets by deployment can be found in Table 5. A comparison of basic statistics, 
based upon the entire 17 months of data before and after the use of the technique 
discussed above, shows very little change in the mean and variances of both time 
series (Table 6). Similarly, the autospectra based on the first year of data alone 
and the 17 month gap filled series, were very similar in the spectral shape with 
little change in the higher frequencies, and more energy contained in the lower 
frequencies of the latter (expected with the longer data set). Based upon these 
results, this technique appears to be have been successful at filling the large gaps 
in the P2 and P3 data sets, while at the same time introducing a minimal amount 
of error into these signals. However, caution must be exercised when examining 
energetics during this period to avoid erroneous conclusions based upon synthetic 
data. 
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Table 5. BASIC DATA STATISTICS FOR MOORINGS P2 AND 
P3: Statistics are based upon six-hourly low pass filtered data for 
actual deployment dates. Units are cm s _1 , °C, and dbar for 
velocity, temperature and, pressure, respectively. Positive 
(negative) values indicate northward (southward) or eastward 
(westward) flows. 



Mooring 


u 




V 


- °v_ 


T 


<3t 


P 


o P 


Deployment 


#1 


































P2 100 m 


-5.43 


11.41 


8.89 


12.58 


9.07 


0.37 


112.14 


0.77 


350 m 


-3.18 


6.78 


7.50 


12.50 


7.05 


0.23 


364.44 


0.78 


500 m 


-3.00 


5.44 


4.09 


12.46 


6.13 


0.20 


519.06 


0.42 


Deployment 


#2 
















P2 100 m 


-6.46 


9.78 


5.45 


13.64 


9.56 


0.23 


131.60 


0.53 


350 m 


-5.36 


6.41 


-3.45 


11.48 


7.13 


0.23 


385.21 


1.26 


500 m 


-3.3 


5.97 


-2.14 


10.72 


5.92 


0.22 


538.09 


1.00 


Deployment 


#3 
















P2 100 m 


-5.05 


10.9 


7.97 


14.28 


9.61 


0.36 


109.26 


0.87 


350 m 


-1.99 


6.73 


5.91 


9.07 


7.06 


0.29 


359.94 


1.17 


500 m 


-0.43 


5.99 


3.20 


9.13 


5.96 


0.18 


510.00 


1.18 


P3 100 m 


-4.65 


10.00 


4.00 


13.34 


9.75 


0.35 


106.60 


1.08 


350 m 


-3.78 


6.76 


2.82 


10.07 


6.88 


0.30 


363.71 


1.56 


500 m 


-3.36 


5.57 


2.52 


8.46 


5.79 


0.22 


513.64 


1.70 


1000m 


-1.57 


3.59 


0.69 


6.70 


3.96 


0.07 


1013.56 


1.69 


Deployment 


#4 
















P2 100 m 


-9.46 


12.56 


14.85 


14.91 


9.52 


0.57 


102.08 


1.70 


350 m 


-3.06 


6.46 


7.62 


9.33 


7.08 


0.30 


350.33 


1.26 


500 m 


-0.52 


6.09 


3.26 


10.27 


5.95 


0.25 


501.64 


0.76 


P3 100 m 


-1.76 


11.24 


2.97 


7.50 


9.56 


6.52 


83.11 


0.95 


350 m 


-1.45 


6.80 


3.42 


5.34 


6.78 


0.27 


346.57 


0.81 


1000 m 


-1.71 


3.71 


0.97 


5.12 


3.94 


0.55 


995.97 


6.38 


Deployment 


#5 
















P2 350 m 


-0.67 


6.55 


4.03 


8.90 


6.95 


0.27 


347.05 


1.62 


P3 100 m 


-4.46 


11.01 


2.62 


11.79 


9.99 


0.61 


80.03 


0.49 


350 m 


-2.84 


6.37 


2.20 


7.68 


6.93 


0.23 


339.20 


1.30 


500 m 


-2.20 


5.81 


1.35 


7.41 


5.88 


0.15 


483.16 


0.70 


1000 m 


-1.36 


4.00 


0.41 


6.34 


3.94 


0.82 


999.53 


1.45 
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Table 6. COMPARISON OF PRE- AND POST-FILL DATA 
STATISTICS AT MOORINGS P2 AND P3: Statistics 
are based upon six-hourly low pass filtered data from May 
1989 to October 1990 for the 100 and 500 m levels and to 
May 1991 for the 350 m level for P2, and from May 1990 to 
May 1991 for P3. Units are cm s' 1 and °C for velocity and 
temperature, respectively, and (cm s _1 , °C) 2 for variances. 
Positive (negative) values indicate northward (southward) or 
eastward (westward) flows. 



Depth 


u 


<5 U 2 


V 


Oy 2 


T 




Pre-Fill 


P2 -100 m 


-6.79 


88.18 


9.68 


165.61 


9.46 


0.19 


P2 - 350 m 


-2.50 


22.18 


4.44 


75.09 


7.04 


0.05 


P2 - 500 m 


-1.61 


13.53 


2.18 


58.00 


5.99 


0.04 


P3 - 350 m 


-2.26 


29.05 


2.71 


25.78 


6.87 


0.05 


Post-Fill 


P2 -100 m 


-6.83 


87.22 


9.51 


163.69 


9.46 


0.19 


P2 - 350 m 


-2.53 


22.37 


4.44 


75.52 


7.04 


0.05 


P2 - 500 m 


-1.60 


13.44 


2.21 


57.68 


5.99 


0.04 


P3 - 350 m 


-2.21 


28.98 


2.65 


25.86 


6.87 


0.05 



3. Time Series Synthesis 

Before one can examine energy conversions and balances through the 
use of the heat equation, it is first necessary to obtain a time series for the 
horizontal gradient of temperature. The horizontal derivatives of temperature 
appear in several terms of the full heat and energy equations, the development of 
which appears later in this text. A time series of VhT can be obtained by using 
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the thermal wind relations as proposed by Niiler and Hall (1988), along with the 
velocity and temperature data from two levels in the array. The resulting time 
series will be for the mid-depth between the two levels used. Therefore, it also 
becomes necessary to synthesize time series of u, v, and T for this same level 
before any further calculations can be made. For moorings P2 and P3 these 
levels will be at 225 dbar and 425 dbar, with an additional level of 750 dbar for 
P3 due to the instrument at 1000 m. All calculations were performed using the 
low pass filtered data sets, because it is within these lower frequencies where the 
thermal wind relation is a useful approximation. 

a. Calculation of Temperatures and Velocities 

There have been several different schemes (see Appendix C) used 
in recent years to correct a time series of temperature for vertical excursions of 
the mooring. In this study, the objective was not to correct the temperature time 
series for mooring motion (a 2 < 2 dbar in most cases) but rather to synthesize a 
time series of temperature for the mid-depth between instrumented levels at 
moorings P2 and P3. Several methods of interpolation (linear, exponential, and 
polynomial functions) were investigated (Appendix C), and while the smallest 
errors were obtained by using a quadratic (cubic) function at P2 (P3) to model T 
vs. P at each site, the errors bars (95% significance level) for each method 
overlapped. This implies that no one method is significantly different from 
another, and therefore the linear approach was selected for use in these 
computations. The equation for linear prediction of temperature can be written 
as 



T(P )new 



T (Ta-TQfr-P) 
P2-P1 



( 1 ) 
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where the subscripts 1 and 2 refer to the upper and lower levels, respectively, 
and P is the pressure level of interest. 

Linear interpolation methods have been used by several 
investigators ( Hall and Bryden, 1985; Hogg , 1986) when correcting velocities for 
mooring motions in deeper water (near 400-500 m); however, they may not be 
suitable for shallower instruments (Hall, 1989b). In regions of intensified 
surface flow, the vertical structure of density and velocity are often modeled 
with an exponential dependence (Gill et al., 1974; Hall 1989b, 1991). The 
velocity calculation scheme used by Hall (1989b, 1991) required the use of 
velocity and pressure records for two levels, and was either a linear or 
exponential interpolation depending upon the signs of velocity at each level. 
When the components of velocity (ui and U 2 ) at both levels were of the same sign 
and not equal, the mid-depth velocity was calculated assuming an exponential 
dependence between the two levels, whereas, when the components of velocity 
(ui and U 2 ) at both levels were of opposite sign or equal in value, the mid-depth 
velocity was calculated assuming a linear dependence between the two levels. 

A problem with this particular scheme is that the results depend 
upon the choice of (x,y) coordinate axes, such that changing the definition of the 
axes changes the type of interpolation used. Additionally, by alternating velocity 
calculations between exponential and linear, it is possible to introduce errors into 
the energy calculations, through false fluctuations which may result from this 
scheme. Because the distances over which we are interpolating are relatively 
small compared to some of these other studies, and the fact that a linear approach 
was found to be satisfactory in the temperature calculations a linear approach 
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will also be used here for the sake of consistency. As with the temperature 
prediction, the velocity prediction equation is formulated as 



u(P) new — Ui - 



(u 2 -ui)(Pi - P) 
P 2 -P 1 



( 2 ) 



where once again the subscripts 1 and 2 refer to the upper and lower levels, 
respectively, and P is the pressure level of interest. Sample time series of actual 
and synthesized data for mooring P3 are shown in Figures 3 through 5, and show 
that this method is acceptable and does not appear to have introduced any 
erroneous signals into the data. 

b. Computation of Horizontal Temperature Gradients 

To compute VhT, the approach of Niiler and Hall (1988) and Hall 
(1991), which utilize the thermal wind relations, has been used. By assuming a 
linearized equation of state of the form p = po(l-otT), the thermal wind relations 
can be solved for horizontal temperature gradients at mid-depth levels as 



where / is the Coriolis parameter (8.5699 x lO 5 s 1 at 36° 20' N), g is the 
gravitational acceleration (980 cm s* 1 ), and a is an equivalent compressibility 
coefficient (j Dewar and Bane, 1985). Vertical velocity shear is determined from 



3T _ J_ 0v 0T _ f Ou 

Ox dz dy dz 



(3) 
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the upper and lower level time series. The value of a is defined as 

a = _ J_dp = _( fx^p] + fx?p]ds + (x^P.)dP ) 

Po dT \\Po dTJ \Po 3S j dT VPo 3P J dT ) (4) 

(a) (b) (c) 

where (a) is the thermal expansibility coefficient, (b) is the salinity contraction 
coefficient, and (c) is the isothermal compressibility coefficient. These values 
have been obtained from the UNESCO International Oceanographic Tables 
(UNESCO, 1987). Based upon the magnitudes of these coefficients, the 
isothermal compressibility term, (c), which is two orders of magnitude smaller 
than (a) and (b), has been neglected in these calculations. The mean value for T, 
S, and dS/dT, at each mooring and synthesized depth, were determined from an 
average T-S profile based upon the CTD data (28 stations) discussed earlier. 
Using this information with the interpolation procedure described in the 
UNESCO tables, a value of a was obtained for each synthesized depth at 
moorings P2 and P3 (Table 7). Once the values of a were obtained, VhT could 
then be calculated through equation (3). The temporal variability of a is small 
and had little effect in the energy calculations. 

The validity of the thermal wind relations in this flow regime can 
be determined by examining the local Rossby number in conjunction with the 
scales of motion expected in this region. The Rossby number, which is the ratio 
of the non-linear terms in the equations of motion to the Coriolis term, is defined 
as 
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( 5 ) 



j£ 

Ro # = -L- = -LL 

/oU f 0 L 

where U, L, and f Q represent the characteristic velocity, length scale of motion, 
and Coriolis parameter, respectively. Table 8 contains values of the local Rossby 
number based upon selected velocity - length scale combinations. Based upon the 
time series data from moorings P2 and P3, peak velocities in the low-pass 
filtered data rarely exceed 40 cm s- 1 , and typical length scales for subtidal 
motions are greater than 50 km, as cited earlier. Therefore, it appears that we 
can expect the local Rossby number to be on the order of 0.1 or less. This is 
consistent with the observations of Pierce et al., (1991), who found the local 
Rossby number to be approximately 0.2 across a CTZ jet, and concluded that the 
regime could still be approximated by QG dynamics. 

At 425 dbar, the values of a at each mooring appear comparable to 
each other; however, this is not the case at the 225 dbar level, where the 
difference becomes larger. The reason for this may come from the fact that the 
value of dS/dT, based upon CTD data, is different at each mooring location. At 
P2 the mean values of T and S are 8.180 °C and 34.077 psu, respectively, while 
at P3 the values were 7.997 °C and 34.075, respectively. The local gradients of 
dS/dT, computed over 20 m, were -0.154 and -0.082, for P2 and P3, 
respectively. If we remove the contribution of these gradients from term (b) in 
equation 2-4 (Table 7), we obtain the nearly identical values of 7.623 x 10' 4 and 
7.622 x 10' 4 , for P2 and P3, respectively. This indicates that the reduced local 
gradient within the mean T-S curve at P3 was responsible for the observed 
difference in the values of a. 
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Table 7. NUMERICAL VALUES OF p and a AT MOORINGS 
P2 AND P3: Computed from UNESCO International 
Oceanographic Tables using an average T-S profile based upon 
CTD data collected during the mooring deployments. 
Components (a) and (b) are the thermal expansibility 
and salinity contraction coefficients, respectively, from 
equation 4. 



Mooring 


Depth 

(dbar) 


p (g cm* 3 ) 


(a)x 1(H 


(b)x 1(H 


a (°C-i)x 1(H 


P2 


225 


1.026524 


-1.506 


-1.174 


2.680 


P2 


425 


1.026850 


-1.377 


-0.443 


1.820 


P3 


225 


1.026547 


-1.485 


-0.626 


2.111 


P3 


425 


1.026863 


-1.356 


-0.527 


1.883 


P3 


750 


1.027206 


-1.255 


-0.803 


2.058 



Table 8. LOCAL ROSSBY NUMBER AS A FUNCTION OF 

LENGTH SCALE AND VELOCITY: The values of the 
local Rossby number (Ro) have been computed in accordance 
with equation 5. Based upon the typical scales of motion 
expected off Point Sur in the CCS and the observed current 
velocities, the value of the local Ro should lie within the 
region to the right of the double line. 



Velocity 
(cm s- 1 ) 


Lengt 


h Scale (km) 


10 


20 


30 


40 


50 


100 


200 


10 


0.10 


0.05 


0.033 


0.025 


0.02 


0.01 


0.005 


20 


0.2 


0.1 


0.067 


0.05 


0.04 


0.02 


0.01 


30 


0.3 


0.15 


0.1 


0.075 


0.06 


0.03 


0.015 


40 


0.4 


0.2 


0.133 


0.10 


0.08 


0.04 


0.02 


50 


0.5 


0.25 


0.167 


0.125 


0.10 


0.05 


0.025 
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Figure 3. U- Component of Velocity at Mooring P3: Data has 
been low pass filtered and decimated to six hourly points. 

From top to bottom are the 100, 225, 350, 425, 500, 750, and 
1000 dbar levels. Data at the 225, 425 and 750 dbar levels have 
been synthesized using data from bracketing levels as described 
in the text. Units are cm s* 1 . 
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Figure 4. V- Component of Velocity at Mooring P3: Same as 
in Figure 3, but for the v-component of velocity. Units are 
cm s' 1 . 
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Figure 5. Temperature at Mooring P3: Data has been low pass 
filtered and decimated to six hourly points . From top to 
bottom are the 100, 225, 350, 425, 500, 750, and 1000 dbar 
levels. Data at the 225, 425 and 750 dbar levels have been 
synthesized using data from bracketing levels as described in 
the text. Units are °C. 
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III. CHARACTERISTICS OF THE POINT SUR DATA SETS 



A. GENERAL OBSERVATIONS 

At mooring P2, the low-frequency current fluctuations between 100 and 500 
m exhibited strong vertical coherence and were characterized by episodes of both 
strong poleward and equatorward flow (Appendix A). The fluctuations were 
distinctly wavelike in appearance, and had typical peak amplitudes of ± 40.0 cm 
s* 1 , a period of approximately 215 days, and were not phase locked with the 
seasons ( Ramp et al., 1991b). Strong poleward flow was observed between 
November 1989 and February 1990, and also between June and September 1990. 
This signal was not apparent in the local wind stress and coastal sea level {Ramp 
et al., 1991b). Currents were, in the mean (based upon all the data), northward 
and westward (north-east coordinate system), and decreased in magnitude with 
depth (Table 6). Similar results were observed using data from individual 
deployments (Table 5), with the exception of the second deployment (August to 
December 1989), where strong shear existed between 100 and 350 m depth 
(Table 5, Appendix A), and the mean flow was to the south and west below 350 
m. While the direction of flow ,was generally the same between deployments, the 
individual statistics did not remain constant indicating the presence of very low 
frequency variability at this location. 

A well known property of autocorrelation functions is that any strong 
periodicities in a given time series will also appear in its autocorrelation function 
( Bendat and Pier sol, 1986). Using this property, it is possible to determine the 
dominant signals within these time series by estimating them as four times the 
first zero crossing of the respective autocorrelation function. The 
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autocorrelation functions for mooring P2 (Figure 6), based upon the same time 
periods listed in Table 6, reveal a very low frequency signal in the alongshore 
currents at all depths, with a period of approximately 230 days at 100 m and 350 
m, and 207 days at 500 m, which is similar to the value of 215 days cited earlier. 
This signal is not present in the autocorrelation functions for the cross-shore 
components which appear to be dominated by fluctuations shorter than 80 days. 
It is interesting to note that the temperature fluctuations at 100 m depth appear to 
be dominated by this 230 day signal, as indicated by the dotted line in Figure 6, 
which closely resembles the autocorrelation function for the alongshore 
component. However, below 100 m depth, the autocorrelation function for 
temperature more closely resembles the cross-shore component of flow, which is 
dominated by shorter period fluctuations. The presence of this very long period 
oscillation will influence the determination of the mean flow at this location, and 
may indicate that future studies will require more than 17 months of data to 
obtain a better estimate of the temporal mean. 

Farther offshore, at mooring P3, flow was also highly coherent over the 
instrumented range of the 100 - 1000 m depths, but was less energetic than flow 
inshore at P2 by a factor of two (Ramp et al., 1991b). On average, the flow was 
poleward with speeds of about 11 cm s* 1 (Table 5); however, the records did 
contain two intense bursts of equatorward flow (velocity greater than 40 cm s* 1 ), 
which lasted approximately 30 to 45 days. These bursts were highly coherent in 
the vertical, and appeared between March and April during both years 
(Appendix A). Whether or not their occurrence during this particular time 
period has any significance or was purely coincidental lies outside the scope of 
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Figure 6. Autocorrelation Functions for Mooring P2: Top, middle 
and bottom panels are the 100, 350 and 500 m levels, 
respectively. Solid (dashed) curve represents u (v) component, 
while the dotted curve is for temperature. 
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this study and will not be addressed here. It is interesting to note that while these 
bursts stand out as the dominant features in the velocity records at P3, they did 
not occur at mooring P2, located only 25 km away. 

The autocorrelation functions for the 100 and 350 m depths at P3 (Figure 7), 
based upon the time period listed in Table 6, reveal that while shorter periods do 
exist, the time series contain a dominant long period fluctuation between 
approximately 120 and 160 days duration. At both levels, all three 
autocorrelation functions appear very similar in structure. 

B. SPECTRAL ENERGY ANALYSIS 

To estimate the frequency distribution of energy present within these records 
at subtidal frequencies, variance-conserving autospectra were employed. 
Initially, the longest time series at each mooring location, namely P2 - 350 m 
depth and P3 - 100 m depth, were analyzed to identify areas of significant energy 
intensification. The time series at the 350 m level at P2 is the longest continuous 
record of this study, and runs from May 1989 through May 1991. At P3, the 
records were considerably shorter, with the longest being at the 100 m level 
between May 1990 and May 1991 (corresponding with data periods used in Table 
6 ). 

The spectra at these two locations were computed using 5 overlapping 
fundamental piece lengths of 243 days at P2 and 3 of 180.5 days at P3. These 
particular piece lengths were chosen to attain the highest possible confidence, 
while at the same time allowing for calculation of spectral energy out to periods 
of 180.5 to 243 days. The records themselves are too short to establish the 
secular variability off Point Sur, and the use of a longer piece length would only 
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Figure 7. Autocorrelation Functions for Mooring P3: Upper 
(lower) panel is the 100 (350) m level. Solid (dashed) curve 
represents u (v) component, while the dotted curve is for 
temperature. 
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act to decrease the confidence in the resulting spectra below significant levels. 
Each data segment was demeaned and detrended, using least squares, prior to 
application of the Fast Fourier Transform (FFT). A Chi-square distribution was 
used to compute 95% confidence limits. In a variance-conserving presentation 
these confidence limits vary with both spectral amplitude and frequency making 
them difficult to display. 

At mooring P2 (350 m), there is a strong peak in the energy spectrum of the 
u-component at the 48 day period (Figure 8), while in the more energetic v- 
component (Figure 9) the peak lies at a slightly shorter period (35 days). Peaks 
at the 27 and 48 day periods occur in the temperature spectrum (Figure 10), in 
general agreement with the velocity spectrum. At longer periods (longer than 80 
- 120 days), the spectra of the u-component and temperature indicate a 
significant reduction in energy. However, the same is not true for the v- 
component, which contains energy from the lower frequency signal which 
appeared in the autocorrelation function (Figure 6). 

Because of fewer data segments (only three segments in one year of data), 
the spectra at P3 (100 m) will have inherently less confidence associated with 
them; however, they do reveal that more energy is contained in the u-component 
(Figure 11) than in the v-component (Figure 12). Peaks in the u-component 
found at the 25 and 90 day periods, contain more energy than peaks in the v- 
component at 90 day periods, as well as at shorter periods. In both the v- 
component and temperature spectra (Figure 13), more energy is contained at the 
shorter periods ( < 20 days), and may be associated with the shallower depth of 
this instrument. 
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Because of instrument failures, most records at P3 were six months or less in 
duration, and only the time series described above provided continuous data for 
at least one year. Therefore, the study of energetics off Point Sur will be limited 
to using the longer time series of P2, which cover a 17-month period between 
May 1989 and October 1990, and the 100 and 350 m records at P3 between May 
1990 and May 1991. 
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Figure 8. Spectral Energy of U-component of Velocity at 350 m 
at Mooring P2: Auto-spectrum of low pass filtered data. 

Plot is based upon two years of data, and has been computed 
using five overlapping pieces as described in text. The 95% 
confidence factors are 0.49 and 3.08, which when multiplied by 
spectral estimates yield the 95% confidence intervals. 
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Figure 9. Spectral Energy of V-component of Velocity at 350 m 
at Mooring P2: Auto-spectrum of low pass filtered data. 

Plot is based upon two years of data, and has been computed 
using five overlapping pieces as described in text. The 95% 
confidence factors are 0.49 and 3.08, which when multiplied by 
spectral estimates yield the 95% confidence intervals. 
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Figure 10. Spectral Energy of Temperature at 350 m at 

Mooring P2: Auto-spectrum of low pass filtered data. 

Plot is based upon two years of data, and has been computed 
using five overlapping pieces as described in text. The 95% 
confidence factors are 0.49 and 3.08, which when multiplied by 
spectral estimates yield the 95% confidence intervals. 
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Figure 11. Spectral Energy of U-component of Velocity at 100 m 
at Mooring P3: Auto-spectrum of low pass Filtered data. 

Plot is based upon one year of data, and has been computed 
using three overlapping pieces as described in text. The 95% 
confidence factors are 0.42 and 4.85, which when multiplied by 
spectral estimates yield the 95% confidence intervals. 
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Figure 12. Spectral Energy of V-component of Velocity at 100 m 
at Mooring P3: Auto-spectrum of low pass filtered data. 

Plot is based upon one year of data, and has been computed 
using three overlapping pieces as described in text. The 95% 
confidence factors are 0.42 and 4.85, which when multiplied by 
spectral estimates yield the 95% confidence intervals. 
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Figure 13. Spectral Energy of Temperature at 100 m at 

Mooring P3: Auto-spectrum of low pass filtered data. 

Plot is based upon one year of data, and has been computed 
using three overlapping pieces as described in text. The 95% 
confidence factors are 0.42 and 4.85, which when multiplied by 
spectral estimates yield the 95% confidence intervals. 
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IV. A CONCEPTUAL MODEL FOR INVESTIGATING 
ENERGETIC INTERACTIONS MEASURED AT A SINGLE 
CURRENT METER MOORING 



In this study, a method of analysis similar to Niiler and Hall (1988) and Hall 
(1991) will be used to investigate the energetic interactions occurring within the 
flow field along the continental slope off Point Sur, California. The eddy 
potential energy equation is used to study these interactions. The basic equations 
and some discussion of their derivation are contained in the following sections; 
however, the complete derivations, along with fully expanded terms, have been 
deferred to Appendix D. 

A. THE THERMAL ENERGY EQUATION AND VERTICAL 
VELOCITIES 

To begin, we consider the basic temperature equation which can be written 
as 



3T dT 0T D n 
— + u — — + v — — + w 0 Z = 0 

dt dx dy 



( 6 ) 



where 0z represents the actual potential temperature gradient at a given level 
(i.e., 225 or 425 dbar in this study). When considering temperature fluctuations 
in the vertical and their role in energy conversions, we must use potential 
temperature, (0), instead of temperature to account for the effects of adiabatic 
compressibility. Using the rules of Reynolds averaging, an equation for the 
mean temperature can be obtained by taking the time average of (6): 
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( 7 ) 




ar -3t 

— — + V — — 



9x dy 



+ v' + W 0 Z + w'0'z = 0 

a y 



- or - 



v h VrT + vh' • V hT" + w 0 Z + w'0' z = 0, 



where the overbar represents the time average of a quantity, and 9T/0t is zero. 
By definition, the time average of a product, AB, is equal to AB + a'b', where 
A = A + a' and B = B + b'. Therefore, the equation for the mean temperature 
includes interactions between the mean and eddy fields, where the eddy fields can 
encompass a wide range of frequencies. By subtracting (7) from (6), we obtain 
the equation for the fluctuating temperature field 



where primes denote fluctuating quantities, the subscript H refers to horizontal 
components and derivatives, and all vertical components associated with 
nonsubscripted terms involve 0' in lieu of T\ Here again the nonlinear advective 
portion may contain contributions from a wide range of frequencies. 

As pointed out by Niiler and Hall (1988), the definition of mean and eddy 
fields in these equations must be considered carefully. The calculation of the 
mean is restricted to the length of the time series, which in this case is 17 
months, yet as seen in the v-component of flow (Figures 6 and 9), energy does 
exist at very low frequencies (longer than 200 days). These very low 



^ + v • VT' + v' H • V h T + w’ 0 Z - v’ • VT' = 0 



( 8 ) 
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frequencies will affect the estimate of the mean and will require careful 
consideration. 

The terms necessary to evaluate these equations at the 225 dbar and 425 dbar 
levels are u, v, w, T, VhT, and a. Details on the procedures used to calculate 
these terms were previously discussed in Chapter II, with the exception of 
vertical velocity, w. Vertical velocity is important in the conversion of potential 
energy to kinetic energy, and vice-versa, through vertical eddy heat fluxes and 
can be calculated at each depth from equation (6) as 



w = - 




3T ^ 0T 

9x v a y 




(9) 



where the local time change of temperature at these depths has been computed 
using a centered finite-difference scheme. This method has been used by Hall 
(1991) to compute vertical velocity in the Kuroshio Extension (bottom depth in 
excess of 5000 m). In the present study, time series of vertical velocity were 
computed using both 0 Z and 0 Z at each mooring (see Figures D1-D3 in Appendix 
D). Based upon these results, it appears that either method could be used to 
obtain vertical velocities at these locations; however, as discussed in the next 
section, terms involving vertical advection of eddy energy (which appear as a 
result of the time variability of 0 Z ) were of comparable magnitude to the 
horizontal advective terms. Therefore, the time series of the actual vertical 
temperature profile, 0 Z , have been used to compute the vertical velocities used 
throughout this study. 
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B. THE EDDY POTENTIAL ENERGY EQUATION 

The equation for eddy potential energy (EPE) can be found by multiplying 
the fluctuating heat equation (8) by pog«T’/0 z . This multiplicative factor is the 
same as that used by Hall (1991), except that it has been multiplied by a mean 
density, po, in order to obtain units of ergs cm-3 s- 1 . This difference comes from 
the fact that a in her formulation had units of gm cm' 3 °C‘ 1 , whereas in this 
study a has units of °C' 1 . Values of po (mean density) have been determined at 
each level based upon CTD data, whereas the value of 0 Z comes directly from the 
time series of 0 Z (calculated directly from time series of 0 above and below) at 
each depth. Defining eddy potential energy as EPE = ^Pog<*T' 2 0 Z , after Hall 

(1991), this equation can be written as 

+ V • VEPE + v' • VEPE + 

dt 

(1) (2) (3) 

vVT' • V h T + pogawT = p °g_ aT ' 7 • VT' 

0 2 0z (10) 

(4) (5) (6) 

where once again overbars and primes denote time mean and fluctuating 
quantities, respectively, and the subscript H refers to horizontal components and 
derivatives. Briefly, the first three terms represent: (1) local growth or decay of 
eddy potential energy; (2) advection of eddy potential energy by the mean flow; 
and (3) self-advection of eddy potential energy by the fluctuating eddies. 

Term (4) represents a conversion of energy between eddies and the mean 
baroclinic flow. Hall (1986a, 1991) interprets downgradient eddy heat fluxes 
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(term 4 above) as a signature of baroclinic instability processes, while 
downgradient eddy momentum fluxes (terms which appear in the kinetic energy 
equation) are associated with barotropic instability. Similar interpretations have 
been made by Dewar and Bane (1985, 1989b), where the eddy potential energy 
equation was derived from the density equation instead of the temperature 
equation. Dewar and Bane (1989b) have also interpreted these processes in 
terms of the direction of eddy momentum and heat fluxes relative to the mean 
momentum and heat gradients, and found that the release of energy to the eddy 
field will occur if the eddy fluxes are downgradient (-v'hT'-VhT > 0). Term 
(5) appears with opposite sign in the eddy kinetic energy equation and represents 
exchanges between eddy potential and eddy kinetic energy. In the ocean, with z 
taken to be positive upwards, a downward heat flux (-w'T > 0) represents a 
conversion of EKE to EPE. This situation occurs when either cold water is 
moved upward or warm water is moved downward, so that isothermal surfaces 
are displaced away from their mean positions. Conversely, an upward heat flux 
(-w'T' < 0) represents a conversion of EPE to EKE, and occurs when either cold 
water is moved downward or warm water is moved upward, so that isothermal 
surfaces are brought closer to their mean positions. Term (6) is a residual term 
based upon the mean eddy advection of the disturbance temperature, a constant, 
and T (the only non-constant quantity in the term) and has no clear physical 
interpretation. It was found to be small in this study and has subsequently been 
neglected. 

It should be noted that the terms involving the spatial variability of the mean 
potential temperature gradients, 0 Z , have been omitted from equation (10). This 
omission is necessary since we only have data from single moorings, which do 
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not allow for determination of this variability. Terms involving spatial 
derivatives of 0 Z come from rewriting the advective term in the eddy potential 
energy equation as follows 

• v(poga^) = v-Vepe - pogavl^-- v(=U 
0 Z ' 2 / 2 \ej (ii) 

(a) (b) (c) 

where (a) is the direct result of multiplying equation (8) by PogOcT'/0 z> an( j (5) 
and (c) come from using the relation V*(ab) = a»Vb + b»Va. As mentioned 
above, we do not have the data set to calculate the actual spatial variability within 
0 Z ; however, as shown in Appendix D, it is possible to utilize data from a 6- 
month period common to both P2 and P3 to examine the effect of neglecting 
these terms. Using data during this common time period, as outlined in 
Appendix D, the ratio (x-component only, since P2 and P3 are east-west of each 
other) of term (c) to term (a) in (11) was 0.0113, which is much less than 1, 
indicating that term (a) can be approximated by term (b) alone without 
introducing any significant error into subsequent energy calculations. This result 
is similar to Hall (1991), who found these were typically an order of magnitude 
smaller than other terms in the energy equation, and therefore also neglected 
them from subsequent analyses. 

While the spatial variability of the mean vertical potential temperature 
gradient, 0 Z , has been neglected in equation (10), it is important to note that the 
time variability of the actual gradient (0 Z ), which appears in both the mean and 
eddy advection terms of (10), has not. Comparison of horizontal advection of 
disturbance temperature to vertical advection of disturbance temperature 
(Figures D4-D6 in Appendix D) in all three time series used in this energy study 
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reveals that, in general, the horizontal component of advection is larger than the 
vertical component. However, there are periods where the magnitudes are quite 
similar, and periods do exist where the vertical advective term is the larger term 
(Mooring P3, Figure D6), and may contribute significantly to the energy balance 
of (10). Therefore, it is not appropriate to neglect the temporal variability of 0 Z , 
a result also found by Hall (1991), who also retained these terms in her analyses. 

To obtain the mean values of the terms in the eddy potential energy equation, 
we simply take the time average of (10), which becomes 

+ v • VEPE + v' • VEPE + 

Ot 

(1) (2) (3) 

v'hT' • V h T + poga^T = 0 

e z (12) 

(4) (5) 

where terms 1-5 have the same meaning as in equation (10), except now they 
represent the time mean values. This is also equation (3.1) in Hall (1991). 
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V. RESULTS 



Using techniques outlined in chapters II and IV, time series of u, v, w, T, 
3T/0X, dT/dy, and 9 Z , along with time series of the 5 terms in the eddy potential 
energy equation (eq. 10, see Appendix E) were computed at mooring P2 for the 
225 and 425 m levels (17 months of data) and at mooring P3 for the 225 m level 
(one year of data). The time mean values for each of the five terms in the EPE 
equation, computed in accordance with equation (12), are therefore restricted to 
these record lengths. Prior to discussing results of these calculations, it is 
important to mention that low frequency signals, including seasonal flow 
reversals which are generally observed within the CCS ( Chelton , 1984; Lynn and 
Simpson , 1987), observed variability in the strength and location of the CUC off 
Point Sur ( Tisch et al., 1992) along with the existence of mesoscale eddies, make 
the determination of the temporal means difficult and can act to widen the 
respective error bars on terms in the energy equation. Such low frequency 
variability does exist in these time series as seen in Appendix A and as illustrated 
by the respective autocorrelation functions and autospectra of chapter II. As a 
result we find in some cases, the error bars for the terms in the eddy potential 
energy balance were quite large and did include zero, thus implying the mean 
values are not statistically different from zero. However, in cases where the 
error bars barely include zero, we may be confident that the sign of the term is 
correct and therefore discuss its particular role in the overall eddy potential 
energy balance in this region. 

Error bars shown and discussed throughout this chapter refer to the standard 
error of the mean. While these error bars provide some measure of confidence 



62 



in the calculated means, they do not represent the 95% significance level, which 
would be far too restrictive in this type of study where relatively weak mean 
flows exist in the presence of much stronger mesoscale variability. For the three 
components of velocity and those energy terms which do not include a product 
with a mean value (i.e., all except terms 2 and 4 in (12)), the standard error of 
the mean was calculated in the usual way, dividing the standard deviation (a) by 
the appropriate number of degrees of freedom (u). For each time series a 
unique number of degrees of freedom was determined by dividing the total 
record length by the value of the first zero crossing of the associated 
autocorrelation function, after Ramp (1989). The standard errors of the 
components of velocity (u, v, and w), VEPE (three components), horizontal eddy 
heat flux (u'T* and v'T), and mean horizontal temperature gradient (VhT) were 
computed in this way. For terms 2 and 4 in equation (12), which involve the 
product of two mean quantities, the error propagation technique of Brooks and 
Niiler (1977) was used. For example, the standard error of the x-components is 
computed as 



O- 

u 



dEPE = <J U C T'T'x I + | U | O c x"p x + G u G c jx x 



dx 



and ' _ 

c^.cfcr ) = c^f c| T x | + 1 u'T’ | o c j x + <VF o c Z 

dx 

where c is a multiplicative constant equivalent to apog0 z 1 that has a unique value 
for each level at each mooring, and the overbars represent a time-mean value. 
In other words, the standard error for each component is the sum of the standard 
error of one quantity multiplied by the absolute value of the mean of the other 
(and vice-versa) plus the product of their standard errors. By using absolute 
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values, we eliminate possible cancellation of terms and obtain the maximum 
possible error from this technique. The net standard error for these terms can 
then be calculated as the sum of squares of the above standard errors as 



r7rr>r V O^v-comn + t^v-cnmn 



V-VEPE 



x-comp T ^ y-comp ^ w z-comp 



and 



0 vhTVh(cT) - V 0 2 x -comp + <3 2 y-' 



comp 



where again the overbars represent the time-mean value. 

A. THE TIME-MEAN EDDY POTENTIAL ENERGY BALANCE 

OFF POINT SUR AT MOORINGS P2 AND P3 

At mooring P2, the largest source of EPE resulted from a net advection 
toward this location (convergence of EPE) by the mean flow (Figure 14, Table 
9) balanced by a net conversion of eddy potential energy (EPE) to eddy kinetic 
energy (EKE) at 225 m, through vertical eddy heat fluxes (term 5 in (12)). The 
eddy field, unlike the mean, represented a sink through a net downstream 
advection (divergence of EPE) of EPE. As cited earlier, the CCS is generally 
considered to be baroclinically unstable due to vertical shear within the water 
column, which would imply that the eddy field should gain potential energy from 
the mean field (i.e., weaken the mean temperature gradient). This was found to 
be the case at 225 m depth, where a somewhat weaker source of EPE comes 
from horizontal eddy heat fluxes acting on the mean temperature gradients (term 
4 of (12)), signaling the presence of baroclinic instabilities within the flow field. 
These downgradient heat fluxes produced a net conversion of potential energy 
from mean (MPE) to the eddies (EPE). As a result of these processes, there was 



64 



TABLE 9. MEAN AND STANDARD ERRORS FOR SELECTED 
TERMS IN THE EDDY POTENTIAL ENERGY 
EQUATION AT MOORINGS P2 AND P3: Values for 
mooring P2 are based upon the 1 7 month period between May 
1989 and October 1990 ( 2055 data points), while those for P3 
are based upon the one year period between May 1990 and May 
1991 ( 1448 data points). The multiplicative constant (c) is 

equivalent to otpogG/ 1 and has a unique value for each level at 
each mooring. The sign convention for the five terms in the 
energy equation is the same as described in chapter IV, where 
the local rate of change of EPE is on the Ihs of equation (12) 
while the others have been evaluated as if they were on the rhs 
of (12), such that a +/- represents a gain (source) / loss (sink). 





P2 


225 m 


P2 


425 m 


P3 


225 m 


Term 


P 




P 




P 




u (cm s* 1 ) 


-5.18 


± 1.23 


-2.61 


±0.80 


-2.69 


± 2.29 


v (cms' 1 ) 


7.38 


± 3.44 


3.49 


±2.89 


2.70 


±2.14 


w (cms -1 ) 


-0.0017 


±0.0008 


-0.0020 


± 0.0009 


0.000 


±0.0010 


T (°C) 


8.39 


±0.13 


6.64 


±0.09 


8.19 


±0.10 


uT (°C cm s -1 ) 


-0.17 


±0.32 


-0.14 


±0.18 


-0.77 


±0.44 


vT (“Corns' 1 ) 


-0.37 


± 0.87 


-0.65 


±0.70 


-0.05 


±0.43 


T x (xlO -8 °C cnr 1 ) 


6.39 


±4.73 


7.78 


±5.25 


0.239 


± 3.03 


Ty (xlO -8 °C cm' 1 ) 


4.62 


±2.12 


5.50 


± 1.95 


1.78 


± 1.76 


e z (xlO -5 °C cm' 1 ) 


9.57 


±0.69 


7.27 


±0.49 


11.4 


±0.88 




x 10' 3 


(ergs 


cm' 3 s' 1 ) 








EPE, 


-0.013 


± 0.055 


+0.010 


± 0.026 


0.000 


± 0.083 


v • VEPE 


+0.148 


± 0.275 


-0.040 


±0.112 


-0.015 


±0.107 


7 • VEPE 


-0.088 


±0.114 


-0.104 


± 0.095 


-0.063 


± 0.048 


v’hT’ • VhcT 


+0.080 


± 0.224 


+0.118 


±0.181 


+0.005 


± 0.076 


apogw’T’ 


-0.153 


±0.060 


+0.037 


± 0.032 


+0.074 


± 0.067 
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Mooring P2 - 225 m 




Local Mean Eddy HEHF VEHF 



Figure 14. Eddy Potential Energy Balance at Mooring P2 - 225 m: 
Mean values of the five terms in the EPE equation, which have 
been evaluated as Local = Mean + Eddy + HEHF + VEHF, such 
that a +/- sign represents a gain (source)/loss (sink). Error bars 
represent the standard error of the mean as described in the text 
Units are ergs cm 3 s* 1 . 



no significant growth or decay of EPE at 225 m depth (the local term is less than 
zero; however, not significantly nonzero). This is an encouraging result since, in 
the mean, it is expected that the total amount of EPE at a given location not grow 
or decay (Hall, 1991). A significantly nonzero value for this term would 
indicate that the sampling period was inadequate for this purpose, and calculated 
means are not representative of their true values. 
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To determine which of the individual components are primarily responsible 
for the sign and magnitude of the values listed in Table 9, we can examine vector 
diagrams of v h and ^hEPE^ and v'hT' and ^h c T (Figure 15). The upper panel of 
Figure 15 illustrates the mean flow advection of eddy potential energy, while the 
individual components and associated error bars are shown in Figure 16 for 
comparison. The mean flow at this depth was toward the northwest (325° T) at 
approximately 9.0 cm s* 1 , while the gradient of EPE indicates that more eddy 
potential energy could be found inshore and slightly north of this location. From 
this we see that the x-component of mean flow advection (+0.213 ± 0.140 ergs 
cm* 3 s* 1 ) was responsible for the energy growth observed at the mooring (Table 
9 and Figure 14) while the y-component (-0.079 ± 0.224 ergs cm* 3 s* 1 ), which is 
nearly orthogonal to the gradient, acted as a sink representing a downstream flux 
of energy. A weak contribution to the growth of energy due to mean flow 
advection also came from vertical advection (+0.014 + 0.074 ergs cm* 3 s* 1 ), due 
to downward vertical flow in the mean and a vertical gradient of EPE, which 
pointed toward the surface. For convenience, the factor of 10* 3 has been 
dropped from these values and those appearing throughout the remainder of this 
chapter. The larger error bars associated with the y-component (Figure 16) are 
to be expected in light of the low frequency signals described earlier and are the 
primary source of uncertainty within the overall advection term. Given the 
existence of these very low frequency signals in the meridional component of 
flow and the fact that the error bars associated with the zonal component of 
advection do not cross zero, it is believed that the resulting sign of the mean 
advection term is correct and represents a source of EPE at this depth. 
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Figure 15. Vector Diagrams of vh, v h epe , v'hT' and v hcT a t 

Mooring P2 - 225 m: Top panel contains v h (solid) and 
V h EPE (dashed, multiplied by 10 5 ) with units as shown. Bottom 

panel contains v^T' (solid) and V H cT (dashed, multiplied by 10 3 ) 
with units as shown. 
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Mooring P2 - 225m 




Total 



x-comp. y-comp. z-comp. 



Figure 16. Components of Mean Flow Advection of EPE at 

Mooring P2 - 225 m: Mean values of the three components 
of the mean advection of EPE, along with the total value from 
Figure 14. A +/- sign represents a gain (source)/loss (sink). 
Error bars represent the standard error of the mean as 
described in the text. Units are ergs cm 3 s 1 . 



The sign of individual components of eddy advection of EPE were the same 
as the corresponding components of mean advection; however, the relative 
magnitudes were not. In this case, the y-component was the dominant term 
having a magnitude of -0.210 ± 0.131 ergs cnr 3 s* 1 , a local loss of EPE which 
represents a downstream growth of energy, while the x- and z-components had 
magnitudes of +0.115 ± 0.123 and +0.006 ± 0.032, respectively. Similar to the 
advection by the mean flow, the resulting sign of the complete eddy advection 



69 



term was the same as that of its largest component, which was significantly non- 
zero. Therefore, we see that while the mean flow acts as a source by advecting 
EPE to this location, net advection in the alongshore direction by eddies 
produces an overall downstream growth of energy, and thus represents a sink for 
EPE at the 225 m level at P2. 

As mentioned earlier, downgradient heat fluxes represent a conversion of 
energy from the mean (MPE) to the eddies (EPE) due to baroclinic instabilities 
resulting from a mean shear in the water column. This process can be seen in 
the lower panel of Figure 15, where the mean horizontal temperature gradient, 
and hence VhcT, implies that a greater amount of mean (available) potential 
energy lies to the north and east of the mooring at this depth (indicating that the 
mean state isotherms slope downward in this direction) and the mean eddy heat 
flux is toward the southwest (205° T) with a magnitude of approximately 0.41 °C 
cm S' 1 . Since these two vectors are in opposite quadrants, both the x- and y- 
components of mean eddy heat flux are downgradient and lead to conversions of 
MPE to EPE. The net values in these directions are of comparable magnitude, 
being +0.031 ±0.123 and +0.048 ±0.187 ergs cm' 3 s* 1 , respectively. If the net 
eddy heat flux vector were upgradient (in the direction of increasing MPE), the 
vectors would lie in the same quadrant and the eddies would then augment or 
feed the mean energy state (MPE) at the expense of local EPE by strengthening 
the mean temperature gradient (increasing isothermal slopes). 

At 425 m depth, the net energy balance is somewhat different than at 225 m 
depth. Once again, eddy advection is represented a local sink, implying 
downstream growth of EPE (Table 9, Figure 17); however, it now represents the 
largest loss term in the balance. In fact, the sign of the x- and y- components are 
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the same and negative (-0.087 ± 0.075 and -0.024 ± 0.086 ergs cnr 3 s* 1 , 
respectively), so that the horizontal eddy field as a whole produces the loss. 
Advection of EPE by the mean flow also represents a loss, although much 
weaker than the eddy term (Table 9, Figure 17), and is dominated by the zonal 
component of flow (Figure 18). The mean flow is again towards the northwest 
(323° T) with a magnitude of approximately 4.4 cm s' 1 ; however, the mean 
gradient of EPE at this depth indicates that EPE increases to the west, and not 
inshore, although it is much weaker than at 225 m. From the orientation of these 
vectors, we can see that the meridional component of the mean flow will not 
contribute much to the overall term because it is nearly orthogonal to the 
gradient of EPE, while the zonal component, on the other hand, is almost 
parallel to the gradient and was, in fact, the dominant component. 

The loss of EPE due to advection appears to be balanced by both horizontal 
and vertical eddy heat fluxes. Downgradient heat fluxes (MPE to EPE) at 425 m 
depth, similar to those at 225 m depth, make this term the primary balance to 
advective loss. The net eddy heat flux is now more southward (192° T) and 
stronger (0.67 °C cm S' 1 ) than at 225 m depth while the ^h°T vector is in 
roughly the same direction (Figure 18). Once again these vectors are in opposite 
quadrants; however, a stronger meridional heat flux is responsible for a greater 
contribution to the total term (+0.090 ± 0.164 ergs cm' 3 s' 1 ) than the zonal 
component (+0.027 ± 0.076 ergs cnr 3 s 1 ), which is nearly opposed to v hcT a 
weaker source of EPE comes from a net downward heat flux (EKE to EPE), 
implying that in the vertical, eddy motions are displacing the isotherms 
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